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A  GRAM-SCHMIDT  SPACE-TIME  ADAPTIVE 
CANCELLER  USING  ADAPTIVE  LATTICE  FILTERS 


I.  INTRODUCTION 


An  adaptive  processor  is  a  filter  that  bases  its  own  design  (its  internal  adjustment  settings)  upon 
estimated  statistical  characteristics  of  the  input  and  output  signals  [1-3].  In  particular,  an  adaptive 
spatial  array  processor  weights  the  coherent  output  from  each  sensor  (antenna)  and  adds  them  to  form 
a  receiving  beam.  For  an  adaptive  array  these  weights  may  not  be  constant  but  rather  can  change  as 
a  function  of  the  spatial  properties  of  the  noise  field. 


The  objective  of  a  space-time  adaptive  canceller  (a  special  form  of  an  adaptive  processor)  is  to 
decorrelate  or  statistically  orthogonalize  the  time  samples  of  a  group  of  auxiliary  sensors  from  a  given 
sensor  (called  the  main  channel).  Figure  1  illustrates  the  general  form  of  a  space-time  adaptive  can¬ 
celler.  From  this  figure,  the  auxiliary  channels  are  linearly  weighted  such  that  the  output  noise  power 
residue  of  the  main  channel  is  minimized,  which  is  equivalent  to  statistically  orthogonal izing  the  auxi¬ 
liary  channels  with  respect  to  the  main  channel.  Consider  a  more  detailed  illustration  of  the  space- 
time  adaptive  canceller  given  in  Fig.  2.  The  main  channel  input  is  designated  by  jc0,  a°d  the  auxili¬ 
ary  channels  are  designated  by  x„,  n  =  1,  2,  ...,  N  -  1.  Let  us  set  Naux  equal  to  N  -  1.  External 
signals  are  received  at  Naux  distinct  auxiliary  sensors  and  sampled  in  time  at  equal  intervals  of  T 
seconds,  using  L  time  delay  taps  per  sensor.  The  main  sensor  is  also  sampled  in  time  but  may  be 
delayed  by  some  arbitrary  time  delay  r. 


AUXILIARY  SENSORS 


Fig.  1  —  Generic  space-time  adaptive  canceller 
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OUTPUT 


Fig.  2  —  Space-time  adaptive  canceller 


We  define  the  nth  auxiliary  input  data  row  vector  x„  as 


x„  =  (*„(0,  x„(t  -  D  ,•  •  •  ,  x„(t  -  (L  -  1  )D) 
and  the  optimal  weighting  row  vector  w„ ,  of  these  L  taps  as 

wn  =  (w’l«.  m’2„.  •  '  '  •  wLn)- 

We  also  define  an  LNaux  length  column  vector,  X,  of  all  auxiliary  input  data  as 

X  =  (x,,  x2,  •  •  •  ,  \N  )T  , 


(1) 


(2) 


(3) 


where  T  denotes  the  vector  transpose  operation,  and  an  LNaux  length  optimal  weighting  column  vector 
W  as 

W  =  (w,,  w:.  •  •  •  .  w,v  )T .  (4) 


We  define  /?„  to  be  the  LN,lllx  x  LNaux  input  covariance  matrix  of  the  auxiliary  inputs  and  r  to 
be  the  LNllliX  length  cross-covariance  vector  between  the  mam  and  the  LNxmx  auxiliary  inputs.  More 
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r  =  £{X'jc0},  (6) 

where  E[- )  denotes  the  expected  value. 

It  can  be  shown  [1-3]  that  W  is  the  solution  of  the  following  linear  vector  equation 

Rxx  W  =  r.  (7) 

One  method  of  finding  W  is  to  use  the  Sampled  Matrix  Inversion  (SMI)  algorithm  [4],  which  is  illus¬ 
trated  in  Fig.  3.  For  this  technique,  the  sampled  auxiliary  covariance  matrice  and  cross¬ 
correlation  vector  f  are  formed,  after  which  W  is  found  by  inverting  R^  and  multiplying  by  r.  This 
algorithm  has  been  demonstrated  to  have  a  very  rapid  convergence  rate.  However,  its  disadvantages 
are  its  numerical  complexity  and  possible  instabilities  related  to  taking  the  inverse  at  R ^ . 


Fig.  3  —  SMI  canceller 


In  this  report,  we  present  a  technique  for  implementing  a  space-time  adaptive  canceller  by  using 
Gram-Schmidt  (GS)  and  adaptive  lattice  filter  (LF)  processing  [5-10],  These  techniques  exhibit  excel¬ 
lent  performance  simultaneously  in  arithmetic  efficiency,  stability,  and  convergence  rate  over  other 
adaptive  algorithms.  The  numerical  stability  of  the  algorithm  is  enhanced  because  it  does  not  require 
the  calculation  of  an  inverse  matrix  as  does  the  SMI  algorithm. 

In  addition,  we  take  advantage  of  the  fact  that  the  time-delayed  samples  in  a  given  auxiliary 
channel  are  statistically  stationary.  That  is,  if  we  compute  the  L  x  L  covariance  matrix,  Rx  ,  ,  of 
the  nth  auxiliary  channel  where 
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then  RX  J(  is  a  hermitian  toeplitz  matrix  that  has  the  form 
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If  the  time  samples  are  stationary  in  a  given  channel,  then  adaptive  LFs  [6-10]  can  be  implemented 
resulting  in  a  significant  reduction  in  the  number  of  numerical  operations. 

In  Sections  II  and  III,  GS  cancellation  and  adaptive  LFs  are  reviewed  respectively.  In  Section 
IV,  the  space-time  GS  canceller  with  adaptive  LFs  (GS/ALF)  '.s  described.  In  Section  V,  the  space- 
time  GS  canceller  is  compared  with  the  GS/ALF. 

n.  GS  CANCELLER 

We  briefly  describe  a  GS  canceller.  Note  that  the  open  loop  digital  processing  of  the  input  data 
is  fundamental  to  this  implementation.  Also  the  auxiliary  channels  to  this  processor  consists  of  all  of 
the  time  taps  of  the  auxiliary  antenna  as  seen  in  Fig.  1. 

For  a  GS  canceller  the  optimal  weights  formulated  by  Eq.  (7)  are  not  computed.  The  data  in 
the  input  channels  are  filtered  directly  through  an  orthogonal ization  network  as  is  demonstrated  in  the 
following  discussion.  However,  the  steady  state  output  residue  power  in  the  main  channel  is  the  same 
as  if  the  weights  were  calculated  exactly  by  the  use  of  Eq.  (7)  and  applied  to  the  input  data  set. 

Consider  the  general  A/ -input  open  loop  GS  canceller  structure  as  seen  in  Fig.  4.  Let 
>0’  '  ‘ '  »  Jm-i  represent  the  complex  data  in  the  0th,  1st,  .  .  .  ,  Af  -  1th  channels,  respectively. 

We  call  the  leftmost  input  y0  the  main  channel,  and  we  call  the  remaining  M  -  1  inputs  the  auxiliary 
channels.  The  canceller  operates  so  as  decorrelate  the  auxiliary  inputs  one  at  a  time  from  the  other 
inputs  by  using  the  basic  two-input  GS  canceller  as  shown  in  Fig.  5.  For  example,  as  seen  in  Fig.  4, 
in  the  first  level  of  decomposition,  yM _ ,  is  decorrelated  with  y0,  V|,  •  •  • ,  \’m _ 2 ■  Next,  the  output 
channel  that  results  from  decorrelating  with  yw_2  's  decorrelated  with  the  other  outputs  of  the 
first  level  GSs.  The  decomposition  proceeds  until  a  final  output  channel  is  generated.  If  the  decorre¬ 
lation  weights  in  each  of  the  two-input  GSs  are  computed  from  an  infinite  number  of  input  samples, 
then  this  output  channel  is  totally  decorrelated  with  the  input:  y  t ,  y2.  ' ' '  •  >'m  -  \- 

Let  y„(m)  represent  the  outputs  of  the  two  input  GSs  on  the  m  -  1th  level.  Then  outputs  of  the 
two-input  GSs  at  the  mth  level  are  given  by 


n  =  0,  1,  •  •  •  ,  M  -m  -  1 , 
m  —  I,  2,  •  •  •  ,  M  —  1 . 


GO) 
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Fig.  4  —  GS  canceller 
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Note  that  y„(,)  =  yn .  The  weight  w„(m\  seen  in  Eq.  (10),  is  computed  so  as  to  decorrelate  yjm  + 1  ' 
with  y^lm.  For  K  input  samples  per  channel,  this  weight  is  estimated  as 


w„(m)  = 


E  y$r?mWy!r\k) 


k  =  1 


where  *  denotes  the  complex  conjugate  and  |  •  |  is  the  magnitude.  Here  k  indexes  the  sampled  data. 


For  the  basic  two-input  canceller,  the  input  to  the  right  is  the  local  main  input  unless  otherwise 


noted. 


III.  ADAPTIVE  LATTICE  FILTERS 

A  transversal  filter  is  a  method  of  processing  time-stationary  sequences.  Its  general  form  allows 
flexibility  in  meeting  or  attaining  specified  filter  requirements  and  has  been  used  extensively  in  digital 
filtering  applications.  Figure  6  shows  the  general  form  of  a  transversal  filter.  Here  each  input  is 
time  delayed  in  equal  increments  (normally  at  the  sampling  rate  of  the  sample  and  hold  circuit).  Each 
of  the  L  input  samples  in  time,  u(t),  u(t  —  T),  .  .  .  ,  u(t  -  (L  —  l)T)  are  weighted.  These  weights 
are  chosen  in  the  optimal  sense  so  that  when  all  of  the  tapped  time  delayed  samples  are  weighted  and 
added  to  the  main  channel,  then  the  main  channel  noise  residue  is  minimized.  Note  that  the  main 
channel  may  be  time  delayed  by  r. 


We  set  u  =  («(/),  u(t  —  T),  ...,«(/-  (L  —  l)7'))r .  Then  it  can  be  shown  that  the 
optimal  weights  are  the  solution  of  the  following  vector  equation: 

=  ruv ,  (12) 

where  Ruu  is  the  L  X  L  covariance  matrix  of  the  tapped  time  delays  and  ruv  is  the  cross-covariance 
vector  of  length  L  between  the  main  channel  and  the  tapped  time  delays. 

If  the  samples  of  u  are  statistically  time  stationary,  then  by  definition 

£(w*(f  -  kT)  u(t  -  jT)\  =  rk_j  ,  (13) 

i.e.,  the  cross  covariance  between  any  two  time  delay  taps  is  a  function  of  only  the  relative  time  delay 
between  the  taps.  The  form  of  Ruu  is  then  a  hermitian  toeplitz  matrix,  where  >\~j>  k , 
j  =  1 , 2,  •  •  •  ,  L  are  the  elements  of  this  matrix. 

Instead  of  solving  Eq.  (12)  directly  for  w,  an  alternative  solution  is  to  orthogonalize  the  u(t) 
input  data  and  then  solve  for  the  optimal  weights.  Let  the  orthogonality  transformation  be  the  L  x  L 
matrix  A.  If  we  set  u'  =  A  u.  then  the  optimal  weighting  vector  w'  must  satisfy  the  following  vector 

equation 

R„  u  ^'  =  IV,  •  (14) 


m 


0± 
‘AVivt 


NRL  REPORT  9054 


AUXILIARY 

INPUT 

0 


r 


u(t) 


v(t) 


u(t-T)  O- 


u(t-2T)  O- 


u(t-(L-1)T)  O- 


l_ 


TIME  DELAY 


J 


OUTPUT 


Fig.  6  —  Transversal  filter 

where  Ru  u  is  the  L  X  L  covariance  matrix  of  the  transformed  tapped  delays  and  ru  l.  is  the  cross¬ 
covariance  vector  of  length  L  between  the  main  channel  and  the  orthogonal ized  inputs.  Because  the 
transformed  tap  delays  are  orthogonal,  the  matrix  is  diagonal  and  hence  the  weights  w'  are 
easily  derived. 

Figure  7  shows  an  implementation  of  the  method  that  statistically  orthogonal izes  and  normalizes 
an  L-tap  input  sequence.  The  basic  two-input  canceller  building  block  shown  in  Fig.  4  has  been  used 
to  orthogonalize  an  N -tap  input  sequence;  i.e.. 


E[u  '**«  )  =  0;  k  *  n;  k,  n  =  1,  2,  ••  •,  L  . 

Because  of  its  form,  the  filter  shown  in  Fig.  7  is  called  an  adaptive  lattice  filter  (LF)  [7J. 


(15) 


We  use  the  adaptive  LF  as  a  basic  building  block  for  the  space-time  canceller  to  be  described  in 
the  next  section. 


Wm 


i 
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i 

I 

}  First  the  L  tapped  delays  of  the  xN  _ ,  auxiliary  channel  are  orthogonal ized  by  using  an  adaptive 

LF,  and  the  proper  weightings  are  found  for  these  L  orthogonal  outputs  such  that  these  outputs  are 
decorrelated  (statistically  orthogonalized)  with  respect  to  the  .r0,  jc  ,,  .  .  .  ,  xN  _2  channels.  In  Fig.  8, 
the  operation  occurs  at  level  1  at  each  of  the  A  -  1  adaptive  LFs. 

Note  that  the  resultant  outputs  from  each  of  the  A  -  1  LFs  on  level  1  are  statistically  time  sta¬ 
tionary.  Hence,  where  initially  we  have  A-  1  time-stationary  auxiliary  antenna  channels  to  decorre¬ 
late  with  the  main  channel,  we  have  reduced  the  number  to  A  -  2  time-stationary  auxiliary  antenna 
channels. 

Next,  the  L  tapped  delays  of  the  output  from  the  far-right  adaptive  LF  on  level  !  is  orthogonal¬ 
ized.  The  L  outputs  are  properly  weighted  to  decorrelate  them  from  the  other  A  -  2  outputs  of  the 
first  level  adaptive  LF.  Hence  after  level  2  of  the  decorrelation  process,  we  have  A  -  3  time¬ 
stationary  antenna  channels  to  decorrelate  with  the  main  channel. 

It  is  apparent  that  the  above  operations  can  be  repeated  until  there  are  no  auxiliary  channels  to 
decorrelate  and  the  main  channel  is  completely  decorrelated  from  the  time  taps  of  the  auxiliary  chan¬ 
nels. 

The  orthogonalization  of  the  time-stationary  data  of  the  rightmost  channel  that  occurs  at  each 
level  of  the  GS/ALF  is  common  to  all  of  the  adaptive  LFs  at  a  given  level.  Hence,  this  need  not  be 
duplicated  for  each  adaptive  LF  at  a  given  level.  A  more  numerically  efficient  implementation  of  the 
GS/ALF  is  illustrated  in  Fig.  9  where  this  redundancy  is  taken  advantage  of  in  the  processing  struc¬ 
ture.  Note  for  this  example  L  =  3,  A  =  4. 

V.  COMPARISON 

Note  that  the  GS/ALF  canceller  is  not  an  implementation  of  the  adaptive  weighting  indicated  by 
Eq.  (7).  However,  it  is  a  canceller  which  uses  less  two-input  building  block  cancellers  and  numerical 
operations  than  the  “full  up”  GS.  Simulation  studies  are  now  being  performed  to  quantify  the 
GS/ALFs  cancellation  and  convergence  performance  and  will  be  documented  in  a  future  report. 

We  compare  the  total  number  of  two-input  building  block  cancellers  needed  for  GS/ALF  imple¬ 
mentation  NCcs//(Zjr  with  the  total  number  of  two-input  cancellers  needed  for  a  “full  up”  GS  imple¬ 
mentation  NCG5.  Here  we  assume  that  there  are  A  antenna  channels  (main  and  auxiliaries)  and  L 
taps  per  channel. 

For  the  complete  GS  implementation  of  this  space-time  canceller  we  can  show  that 

NCgs  =  |  [(A  -  1  )L  +  1]  (A  -  1  )L  ~  y(AL)2.  (17) 

For  the  GS/ALF  implementation,  we  can  show 

KCgs/alf  =  (A-  l)[L  (,5A  +2)  — 3]  =  yA2L.  (18) 

Hence 


NCGS/alf  _  _1_ 
NCcs  ”  L  ' 


(19) 


•  '  denotes  main  input  to  canceller 


T 

MAIN  RESIDUE 

Fig.  9  —  Efficient  implementation  of  GS/ALF 


or  the  GS/ALF  canceller  implementation  uses  approximately  ML  fewer  building  block  cancellers  than 
the  “full  up"  GS  canceller  implementation.  For  example,  if  there  are  live  tune  taps,  the  GS  Al  l 
uses  approximately  one-fifth  as  many  building  block  cancellers  as  the  GS  canceller 

It  can  also  be  shown  that  the  numerical  efficiency  (total  number  of  complex  multiplications'  of 
GS/ALF  canceller  to  the  GS  canceller  is  also  inversely  proportional  to  the  number  ot  tune  taps 
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